Excitation spectrum for an inhomogeneously dipole-field-coupled superconducting qubit chain 
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When a chain of A'^ superconducting qubits couples to a coplanar resonator in a cavity, each of the qubits 
experiences a different dipole-field coupling strength due to the waveform of the cavity field. We find that this 
inhomogeneous coupling leads to a dependence of the qubit chain's ladder operators on the qubit-interspacing 
I. Varying the spacing I changes the transition amplitudes between the angular momentum levels. We derive an 
exact diagonalization of the general A^'-qubit Hamiltonian and, through the TV = 4 case, demonstrate how the l- 
dependent operators lead to a denser one-excitation spectrum and a probability redistribution of the eigenstates. 
Moreover, we show that the variation of I between its two limiting values coincides with the crossover between 
Frenkel- and Wannier-type excitons in the superconducting qubit chain. 

PACS numbers: 02.20.Qs, ()3.65.Fd, 42.50.Ct, 85.25.-j 



I. INTRODUCTION 

Superconducting quantum circuits have attracted consider- 
able attention because of their capabilities (i) to demonstrate 
macroscopically the basic interaction of one "atom" and one 
photon in a cavity (e.g., Ullll])- (ii) to serve as aplatfrom for 
testing many quantum optical phenomena (e.g., Ij-uIDj as well 
as (iii) to show its potential as a basis for quantum information 
processing (e.g., Isl-lllll). 

While research on single-qubit interactions are more com- 
mon, many recent articles also studied multi-qubit interac- 
tions. Superconducting circuits with such interactions are also 
known as quantum metamaterials lfl2ll . To be precise, the cir- 
cuit system we consider here consists of a chain of N super- 
conducting two-level qubits coupled to a photon mode in a 
superconducting coplanar waveguide resonator. Compared to 
the single-qubit version, it can manifest even more quantum 
phenomena, including plasma waves ifisll . controllable col- 
lective dressed states \\%, and quantum phase transitions ifTsl - 
[Tsll . It also promises potential for various applications, includ- 
ing quantum simulators ifTgll and quantum memories ll20ll . 

Theoretical studies of multi-qubit interactions with a pho- 
ton often employ the Dicke model 121], where the Pauli op- 
erators are summed and transformed into a bosonic operator 
In this approach, the chain of qubits is treated collectively as 
an atomic ensemble and the excited qubits are collectively re- 
garded as one exciton mode. This theoretical simplification 
proves adequate when (i) the number of excitations in the sys- 
tem is low (in the so-called "one-photon" processes) and (ii) 
the number of qubits is large enough such that the interspac- 
ing Lq between neighboring qubits can be ignored compared 
to the photon wavelength Lp in the resonator (i.e., the qubits 
can be regarded as a continuum). 

However, the question of how excitations arise in supercon- 
ducting metamaterials when these two conditions are not met 
remains unanswered. In a realistic setting for a superconduct- 
ing circuit, the number N of qubits present can range from 
one to, say, 10, but N would not be as large as the number 
of atoms we usually have for an alkaline atomic ensemble in 



an optical microcavity, which is typically greater than 10^. 
Therefore, the Dicke model, which treats N ^ oo, does not 
apply well to the case of multi-qubit superconducting circuits 
with N < 10. 

When a chain of superconducting circuit qubits is arranged 
as a one-dimensional array (i.e., a superconducting qubit 
chain or SQC), each qubit is inhomogeneously coupled to the 
circuit photon mode. In other words, each qubit has a different 
coupling strength to the traversing photon field. This occurs 
naturally since, unlike its optical cavity QED counterpart, the 
photon wavelength Lp is comparable to the qubit interspac- 
ing Lq in a superconducting circuit. The coupling strength 
thus depends on the position of the qubit relative to the pho- 
ton waveform. The effect of the varying coupling strength be- 
comes even more obvious if multi-mode couplings are taken 
into consideration. For example, the qubits on the antinodes 
of the waveform will couple most strongly, whereas those on 
the nodes will not couple. 

The first step to understand and characterize this 
inhomogeneously-coupled system (the aim of this article) is to 
obtain the energy spectrum of the collective excitation mode 
in the chain of qubits and to compare it with that of the Dicke 
model. We find that the inhomogeneity of the couplings in- 
curs a n algebr aic deformation of the Pauli operators of the 
qubits 1I22I - I25I1 . We quantify this deformation through a "de- 
formation factor," which is a function of the relative spacing 
I = 2Lq/Lp, and characterize the amount the inhomogeneous 
system deviates from the homogeneous case. The deformation 
factor modifies the spin operators of the collective qubit chain. 
Consequently, the excitation spectrum will not only be a func- 
tion of the eigenenergy of the photon mode and the qubit level 
spacing, but is also highly related to the deformation factor 
and hence the relative spacing I. 

Note that when atoms are confined to a cavity, the mag- 
netic or laser field that is exerted on them is uniform. The 
strength of the interaction can be uniformly increased or de- 
creased according to the density of the atoms. This macro- 
scopic viewpoint does not differentiate between the identities 
of the atoms. However, for circuit QED, the identities of the 
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qubits are partially differentiated since the qubits can be cat- 
egorized according to the values of their coupling strength to 
the photon mode. This partial differentiation has made un- 
derstanding the inhomogeneous system a many-body physics 
question. 

Our deformation algebraic approach here is a statistical ap- 
proximation method that can be regarded as finding the aver- 
age contribution of the coupling strength given by the SQC 
as a whole. In the end, the characterization (the excitation 
spectrum) of the SQC as an inhomogeneous system is not 
parametrized by the individual qubits, but by the relative spac- 
ing I. In other words, the spacing I is one extra degree of 
freedom peculiar to the inhomogeneous SQC, not seen in a 
homogeneous optical cavity. 

We will first introduce the model and derive the deforma- 
tion factor in Sec.|ITl With the deformation factor, new opera- 
tion rules for the spin angular momentum operators are found 
by solving a difference equation in Sec. |lll] The general en- 
ergy spectrum for n-qubit SQC is given in Sec. UV] We also 
derive in Sec.|IV]a one-excitation spectrum for a 4-qubit SQC 
as a nontrivial case to show the effects of the inhomogene- 
ity. Namely, the energy splittings between the eigenstates of 
the deformed coupling case shrink, while the probability am- 
plitudes of the eigenstates are redistributed such that higher- 
photon occupations are favored. 

In the final Sec.[V] we will consider how the collective ex- 
citations on the SQC would emulate the excitons in atomic 
lattices. In one limit, it becomes a Wannier-type exciton ll26ll . 
where the wave function of the excited level is localized on a 
single atom. In the opposite limit, it emulates a Frenkel-type 
exciton, which has an extended wave function across multiple 
atoms. Changing the degree of freedom I lets the emulated 
exciton undergo a crossover between these two types of exci- 
tons. This crossover depends on a (27V — l)-th order trigono- 
metric equation, whose solution corresponds to the asymptotic 
turning point from the deformed (inhomogeneously-coupled) 
SQC to the undeformed (homogeneously-coupled) SQC. 



II. INHOMOGENEOUS COUPLING MODEL 
A. Inhomogeneous coupling 

For a finite number N of spins in the SQC, the prob- 
lem discussed here is similar to the Tavis-Cummings (TC) 
model ll27l l28ll . where all the spins are grouped into a total 
"large" spin. However, the exactly solvable TC model ap- 
plies only when the coupling is homogeneous and when the 
eigenfrequencies between the qubits and the photon mode are 
equal. When the coupling is inhomogeneous, the large spin 
does not obey the usual commutation relations of the Pauli 
matrices, which the TC-model assumes. 

The new commutation relations of the large spin introduced 
by the inhomogeneity are pertinent to the deformed SU(2) Lie 
algebras. From these algebraic structures, we can establish a 
deformed dipole-field coupling model, of which the TC model 
is a special case. In the following discussion, we consider the 
typical case where the inter-qubit spacing is uniform. The 




Figure 1 : (Color online) (a) Schematic diagram of an SQC with spac- 
ing 1 — 1/3. The upper strip represents a coplanar resonator whose 
potential waveform of wavelength Lp is drawn as a sinusoidal curve, 
while the lower one represents a ground strip. The rectangles be- 
tween the strips are qubits with interspacing Lq. The gray dots on 
the curve indicate the anti-nodes of the waveform, (b) The elonga- 
tion effect on the large spin due to inhomogeneous coupling. The 
ellipsoid shows a case of deformation -Rjvj = 0.5. The unit sphere 
shows the spin under the usual homogeneous coupling. 

coupling strength of each qubit to the photon field can, there- 
fore, be written as a cosine function of a phase factor which 
is determined by the position jl of the j-th qubit, where j is 
the reduced coordinate and I is the relative spacing introduced 
above. 

The situation is illustrated in Fig. [Ha). A chain of qubits is 
sandwiched between the superconducting coplanar resonator 
and a superconducting ground strip. The photon mode pro- 
viding different potential energies on the spins is shown by 
the red sinusoidal curve. 

We use the operators {(Jj,z} to denote Josephson junction 
qubits, and {a,a^} to denote the operators for the single- 
photon mode. With the wave vector being the reciprocal 
of the photon field wavelength on the one-dimensional lat- 
tice, k — 27r/Lp, the dipole-field coupling is of the form 
aj,x{a cos(j7rZ) + h.c). Under the rotating wave approxima- 
tion, the Hamiltonian can be written as {h = c = 1) 

H — u!q (7j^z+woa^a+77 cos(j7rZ) [aj^+a + aj^-a^] 

J=0 3=0 

(1) 

where Wq is the eigenenergy of the spins, ujq the mode fre- 
quency of the photon, and rj the coupling amplitude. 

To diagonalize the Hamiltonian in Eq. ([T]), we introduce the 
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"large spin" operators: the magnetic moment or z-direction 
collective spin operator 



N-l 



(2) 



which is no different from the homogeneous case, and the 
paired raising and lowering operators 



JV-l 

S+ = ^ cos(j7r/) 

N-l 

S- = ^ (Tj,_ cos(j7rZ) 

3=0 



(3) 



(4) 



which have the special sinusoidal dependence on I due to the 
inhomogeneity. The commutator of the paired ladder opera- 
tors no longer equals to 2Sz but has an additional term due to 
the cosine coefficients, i.e., [5*+, S-] — 21]^ with 



N-l 

= S'z + ^ siii[j7r(l 



(5) 



The detailed derivation is shown in Appendix I A II Note that 
in the usual circuit QED system where only one spin is 
placed at midway, the spacing I equals to 2, for which the 
latter term in vanishes. This is the limiting case which 
corresponds to the Wannier type of excitation, where the set 
of spin operators retains the usual structure of an undeformed 
SU(2) algebra. 



B. Deformed algebraic structure 

When the second term of does not vanish, the algebraic 
structure is called deformed ll22l - i25ll . In order to quantify the 
deformation, the commutator of the ladder operators needs 
to be expressed as a function of Sz, i.e., Sz = ,f{Sz). To 
find this function /, we consider an underlying manifold, on 
which there is a local point, say the origin 0, where we define 
a tangent space with the Pauli z-matrices {<Jj.z} being its ba- 
sis vectors, since these matrices are linearly independent. The 
operator Sz, defined above with uniform coefficients, can be 
deemed a vector in this tangent space; the operator is then 
another vector dependent on the parameter / and is a deviation 
or deformation from Sz- Thus, the first-order approximation 
of Sz with respect to Sz is its projection onto the vector Sz- 
That is, since the cosine coefficients are bounded, we can use 
their Hilbert-Schmidt norm 

(Sz, 5z) = tr (S;^z) = N+Y^ cos{2jTTl) (6) 

j=o 



and the Schmidt decomposition 112911 to write Sz = Rn,iSz as 
a deformation of the original z-spin operator where 

1 r sm\(2N ~ 1)t:1]\ 

Rj,^i = \2N + 1+ ^\ , ^ (7) 

4:N I sm(7r/) 



is the deformation factor (Cf. Appendix lA 21 for this deriva- 
tion). The commutator of the ladder operators can now be 
expressed as 



[S+, S^] — 2Rn^iSz, 



(8) 



where Rn.i has a limiting value of one when / — > or / — > oo, 
for which the usual structure used in the TC model is retained. 

Since the deformation factor Rn.i does not affect the com- 
mutation relations between the ladder operators and the z- 
spin, the large-spin operators {<Sz,5'+,5_} form a specific 
deformed algebra ll22[l23ll and not the more general type ll24ll . 
The Casimir operator 



C = S-S+ + h{Sz) 



(9) 



of the algebra, which equals to the undeformed spin momen- 
tum square, S^ ~ S"^ + Sy + Sz, for the homogeneous 
coupling case, is accordingly deformed. Through solving a 
recursive relation (Cf. Appendix |A 31 for details), we find 

h{Sz) = Rn,i{SI + Sz) and hence 



C — S*^ + S*,^ + Ri\[,iS'^, 



-■x ' 



(10) 



which shows that the spin momentum is reduced along the 
z-direction; 



D — + 0„ + ltj\[ lOz 



(11) 



To visualize this reduction of the spin momentum, we can take 
a unit value for the spin moment and let the Casimir operator 
be represented by a unit sphere in 3-dimensional space for the 
undeformed case. With Rn.i < 1, the deformation would be 
an elongation of the unit sphere, along the z-axis, to an ellip- 
soid, while the x- and y-semi-minor axes remain unchanged, 
as shown in Fig.lHb). 



III. OPERATION RULES 

If we consider the unit sphere of Fig.[TJb) as a Bloch sphere 
on which the large spin prescribes its N levels, its elongation 
due to deformation will accordingly modify the transitions be- 
tween the levels. Since the spin-up and the spin-down mo- 
menta do not change, which are still ujqN/2 and —ujqN/2, 
the narrow part of the ellipsoid effectively squeezes the tran- 
sition probabilities. 

More precisely, we consider an arbitrary eigenstate \r,m), 
for which 



Sz \r, m) = m \r, m) 
S'^\r,m) = r(r + l)\r,'m). 



(12) 
(13) 



The ladder operators result in an (r, m) -dependent off- 



(r) 

diagonal matrix element «,„ , i-C-, 

S+\r,m) = a^^^jr, m + 1) 
S- |r,m) = a^'^^i |r,m- 1) 



(14) 
(15) 
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By examining the diagonal elements of the commutator of 
the ladder operators, we find a difference equation (am')^ ~ 
(a^^Li)^ = -2mi?A,_/. With the value a^ll = 0, the equa- 
tion can be solved to give the deformed off-diagonal matrix 
elements or transition probabilities for the ladder operators 

= ^RN,iir-m)ir + m+l). (16) 

See AppendixiBjfor its derivation. 

Geometrically speaking, the deformation process is a 
homeomorphism with a redefined metric g = {1,1, Rnj). 
Since Rn.i < 1, the metric norm is less than unity. The de- 
formation does not affect the level spacings of the magnetic 
moment S^' the number m still takes (2r + 1) values (i.e., the 
ellipsoid is homeomorphic to the sphere). But the transition 
amplitudes to traverse the sphere decrease: if we start with a 
spin-up state \r,r) and finish with a spin-down state |r, — r), 

then all iterations with S- \r,m) ~ Q^mLi \r,m— 1) have 

smaller than the original o;™^ = \/ {r — ■m){r + to + 1) 
(i.e., the ellipsoid is not isometric to the sphere). 

The deformation factor expressed in Eq. d?) is an oscillat- 
ing function of I, where the sine in the numerator determines 
the period of oscillation and the sine in the denominator de- 
termines the period of the envelope. Therefore, the spin an- 
gular momentum of the SQC would be oscillating between 
the unit sphere and the ellipsoid, depending on the qubit spac- 
ing /. The plot of Rn,i in Fig. Eta) for an SQC of = 30 
qubits shows a typical case with envelop of period 1 and local 
minimum of 0.4. The function h{Sz) associated with this de- 
formation factor is a parabola of the magnetic moment 5r. 
For a nontrivial deformation Rn.i < 1, this parabola flat- 
tens and the spin levels become denser The curvature of 
the parabola decreases while its minimum value —R^j/A in- 
creases. As shown in Fig. |2b), the black (gray) arrow indi- 
cates the spin level r = 1/2 for the deformed Rn,i = 0.4 
(undeformed Rn,i = 1) case of SQC. So varying I makes the 
curve h{Sz) oscillate between the boldened curves that cor- 
respond to Rn,i = 0.4 and Rn,i = 1, respectively. We can 
also observe that the level splittings are reduced, reflecting 
the elongated structure of the Bloch sphere in Fig. [Ttb) and 
the modified operation rule of Eq. (fT6T l. 



IV. DEFORMED SPECTRUM 

Equipped with the modified operation rules, we can diago- 
nalize the Hamiltonian in Eq. ([T]i. First, we split the Hamilto- 
nian into two parts 

Ho = uo{S,+a^a), (17) 
Hi = woata + ?7(5+a + 5_at). (18) 

Let u be the number of total excitations and hence the eigen- 
value of Hq. Let n be the number of photons in the system 
such that the SQC magnetic moment is m = li — n. Let i/ be 
the eigenvalue of the interaction part Hi , where uj^ ~ ujQ—ujq. 



(a) 




12 3 

Spacing / 
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Figure 2: (Color online) (a) The deformation factor Rn,i\n=30, as 
an oscillating function of the qubit spacing /. (b) Plot of the function 
h{S~) for various values of -Rjvj. The top and bottom boldened 
curves correspond to -Rjv,! = 1 and 0.4, respectively. The gray- 
shaded double-arrows indicate three spin levels of the undeformed 
SQC, while the black arrows show the r = 1/2 level of the deformed 
case. 



The eigenvector of the Hamiltonian can then be expanded 
as a superposition of different configurations of photon num- 
ber and spin states: 

\u,r) ='^Cn\n;r,u-n) . (19) 

n 

The expansion coefficients c„ satisfy a recursive relation ll28ll : 

Cn+Wn+ la*J'2(„+i) - c«f'« + Cn-iVna^ln = (20) 

where €'„ = (w — ujQn)/!]. 
The solution reads 

L»/2J 

E (21) 

p=0 

'T^n,p can be regarded as a probability amplitude contribution 
to the n-photon state from a set of corresponding qubit chain 
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Table I: Configurations \n;r,m) and probability amplitudes c„ for 
the one-excitation spectrum of an SQC with A'^ = 4 qubits and spac- 
ing I = 2/3. Each o indicates one photon while ^ or ^ denotes the 
spin state of each qubit. Here, v„ — (v — toon)/!^. 



u = 1 
r = 2 



photon 

spin 
config. 



n = 
m = 1 



n = 1 
m = 



n = 2 
m = —1 



n = 3 
m = -2 



nn 
mt 
tm 



tm,un 
mt,;m 



an 



1 



States indexed by p\ 



■E 



n 



+ 1) 



■ {ji---Jfc---iLn/2j ) '^=1 

(22) 

where P„ = n"=o and ( ji . . . . . . jL„/2j ) 

represents an index set of descending order {Vfc < ^ : < 
ji < jfc ^ 2; < ji < n — 2}. We can see from Eq. (l2Tl i that 
the operation rules discussed in the preceding paragraphs have 
made the probability amplitudes deformation-dependent, thus 
/-dependent. This will consequently lead to a redistribution of 
probabilities for different photon states. See Appendix IC 1 I for 
the derivation of these coefficients 'T^n,p- 

The simplest nontrivial example of this deformation effect 
can be seen in Table. HI where we consider the one-excitation 
(u = 1) spectrum of a four-qubit SQC with spacing 1 = 2/3. 
The deformation factor in this case is Rnj = 5/8. For a 
weakly coupled SQC with |a)ol ^ the eigenenergies for the 
four levels are given by (Cf. Appendix IC 2l for the derivation) 



E. 



3^ +1 



5u}q ± 4cZ>o i 



36i? 



,1/2- 



1/2 



(23) 

For a deformation factor Rn.i < 1, the splittings between 
these dressed levels are suppressed. In addition, if we substi- 
tute the value of v into the coefficients c„, we will find that ci 
is greater than that of the undeformed case, C2 (ca) increases 
by a greater proportion than ci (C2), while cq remains equal to 
one. Therefore, the probability distribution shifts toward the 
end that favors states with greater number of photons and less 
degeneracy. 



V. EXCITON CROSSOVER 

If the spacing I is large (Lq ^ Lp), i.e., there are multiple 
photon wavelengths between two neighboring qubits, we can 
consider the excitation that the dipole-field coupling induces 
on a qubit to be localized on that qubit. Consequently, this 
type of excitation emulates the Wannier exciton on an atomic 
lattice. If I is small {L^ <C Lp) with a single-photon wave- 
length extending over all the qubits, we can consider the ex- 
citation on the SQC to be delocalized. This type of excitation 




10 10 10 

Reciprocal of spacing l/l 

Figure 3: Semi-log plot of the deformation factor -Rat,; versus l/l 
over one period, showing the Wannier limit on the left end and the 
Frenkel limit on the right end of the horizontal axis. The number of 
qubits is set to iV = 1000. 



emulates the Frenkel's exciton model for molecular crystals. 
Note that when I tends to either zero or infinity in Eq. (|5]), 
falls back to Sz and the regular SU(2) algebra for the commu- 
tators is obtained. Thus, it is justified that in the large- limit, 
the low-energy excitation becomes bosonic for both Wannier 
and Frenkel excitons. 

To determine when the emulated exciton crosses from 
Wannier- to Frenkel-type, we plot in Fig. [3] the deformation 
factor versus the reciprocal of the spacing over 10~^ < / < 1, 
where we can see the asymmetry between the left edge {l/l ^ 
0) for the Wannier limit and the right edge {l/l —> 00) for the 
Frenkel limit. Setting dRN,i/dl = 0, we obtain the trigono- 
metric equation 



tan[(27V - 1)ttI] = {2N - 1) tan{nl). 



(24) 



Transforming Eq. ( l24l i to C/2Ar-i(cos7r/) = 2N ■ 
T2N-i{cosTrl), where T2N-1 {U2N-1) is the Chebyshev 
polynomial of the first (second) kind, we can observe that 
it is a {2N — l)-th order polynomial equation. Hence, the 
curve has {2N — 1) local extrema in exactly {N — 1) oscil- 
lations from the Wannier end to the asymptotic Frenkel end. 
Between these two limits, the excitation has various degrees 
of deformation and the crossover is continuous. We can re- 
gard the crossover point to be the absolute minimum before 
the deformation factor asymptotically approaches one. This 
point approaches when iV — > 00. For the case illustrated 
in Fig. [3] with N = 1000, a numerical estimation gives the 
crossover at / = 7.16 x 10~^, or a length of 2800 spins per 
photon wavelength. 
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VI. CONCLUSION 



2. Deformation factor 



We have studied the inhomogeneous coupling between a 
SQC and a superconducting coplanar resonator, which leads 
to a set of deformation-dependent operation rules of spin mo- 
mentum. The modified rules correspond to tighter energy 
spacings and a shift of the probability distribution of spin lev- 
els. The inhomogeneous coupling also gives rise to different 
types (Frenkel and Wannier) of collective excitations on the 
SQC and the crossover between these types is determined by 
a polynomial equation of the qubit spacing I. 

FN acknowledges partial support from LPS, NSA, ARO, 
DARPA, AFOSR, NSF grant No. 0726909, JSPSRFBR con- 
tract No. 09-02-92114, MEXT Kakenhi on Quantum Cy- 
bernetics, and Funding Program for Innovative R&D on 
S&T. YXL acknowledges support from the NNSFC grant 
No. 10975080 and 61025022. 



Appendix A: Derivation of the deformation factor 

1. Commutation relation 

The commutator of the ladder operators can be computed 
as follows 

iV-l 

[S'+,S'-] = ^ COs{jTTl) COs{kTTl)[aj^+,<Tk,-] 
j.k=0 
N~l 

= 2 ^ cos^(j7r/)crj,2 

j=o 

N-l 

= Sz + ^ cos(2j7rZ)CTj^2. 

To use a consistent notation, we write the R.H.S. as 22^ and 
extract 5*2 from the second term 

2Y.Z = 25, + ^ [cos(2j7r0 - cos(2j7r)] CTj,^ 

JV-l 

= 25*2 + 2 ^ sin [j7r(l + I)] sin [j7r(l - I)] a-j^^ 

3=0 

which gives Eq. (|5]l. 

We can also check that other commutation relations are pre- 
served 

[5'2,5'i] = ^ [(Tk^z,(Jj,±COs{jTTl)] 

N 

= ±(Tj_± cos(j7rZ) 

= ±s±, 

from which we conclude that the newly defined operators 
{^z, 5+, S*-} form a Poly chronakos-Rocek type of deformed 
SU(2) algebra. 



First, we recognize that each qubit tTj ^ has two orthonor- 
mal basis vectors {\ej) , \gj)} for a fixed relative coordinate 
J and hence {|efc) : |efc)'e {|ej) ,\g,) J = 0, . . . , N - 1}} 
forms an orthonormal basis for the Hilbert space H that spans 
all the qubits on the chain. Then Sz for the original spin angu- 
lar momentum and for that of the inhomogeneous SQC be- 
come operators on this Hilbert space H. Since the sinusoidal 
functions cos(j7rZ) are bounded, we can define the Hilbert- 
Schmidt inner product as 



2JV 



(^z,Sz 



k 

2N N~l 

X] ^^'^l X! cos^ U'^O'^j.z ■ cri,z \ek) 

k j,l=0 
2N 

X (efel cos^([fc/2j7rOcr[fc/2j,2 • cr[fe/2j.2 kfc) 

k 

2N 

Xcos2([fc/2j7rO 



k 
N 



= Y,2 cos" {j^l) 
j 

which equals to Eq. (|6]l. We can then write the approximation 
of E, as a Schmidt projection on 5*2 



{-^z.Sz) 



Sz 



{Sz 5 Sz) 

N+1 + Y.J cos{2jTTl) 



2{N+l) 



Sz 



1 
2 

= RnaSz 



N-l 



^ cos{2jTTl) 



Sz 



where Rn.i denotes the deformation factor Its expression can 
be further simplified to 



1 

47V 
1 

m 



2N + 1- cos{2NttI) + 



2N + 1- cos(2Af7r/) + 



sin(27r/) sin(2iV"7rO 



1 — cos(27rZ) 
cos(7r0 sin(2A^7r/) 



sin(7r^) 



2iV + 1 



sin[(27V - l)7rZ] 



siii(7r^) 



where the first line is derived by comparing the real parts in a 
summation of exponentials. 



3. Casimir operator 



The Casimir operator for the algebra is 

C = S-S+ + h{Sz) 
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where the second term satisfies a recursive relation ll24l] 

h{S,) - h{S, - 1) = 2Rn,iS,. 

This relation leads to a solution composed of Bernoulli poly- 
nomials 

h{S,) = RnAB2{-Sz)~B2) 
— Rn,i {S1 + Sz) 

where B2{—Sz) is the second-order Bernoulli polynomial 
with the operator Sz as variable and B2 is the second 
Bernoulli number. The Casimir operator becomes then 

= —{S-^-S- + S-S+) + Rjy^iS'^ 

which equals to Eq. ( fTOl i and represents a deformed total spin 
operator S^. 

Appendix B: Deriving operation rules 

Assume the eigenstate of the z-spin momentum operator 
Sz to be |r, to), that is, r{r + 1) denotes the total spin number 
and m the magnetic moment, for which Eqs. (fT2li-(fT3]l are 

(r) 

satisfied. Further, assume am to be the coefficients when the 
ladder operators are applied to the state vectors as in Eqs. (fT4] l- 
(flSl l, which is indexed by r and to. 

By applying the vector |r, to) to the commutation relation 
Eq. dH), we find 

{r,m\[S+, S-]\r,m) = (r,m\ S+S- — S-S+\r,m) 

= {c.i:i,r - {a<^jr 

= {r,m\2RN,iSz\r,m) 

= 2i?ArjTO. 

Selecting the second and the fourth line, we arrive at a differ- 
ence equation of m: 

To solve the equation, we list out the iterations until the last 

(r) 

entry where Q!_^-i = since —r is the minimum value m 
can take as the magnetic moment 

{a'il^f ~ {a^::l2? = -2i?Ar.i("i-l) 

{a^llf {a^:l_,f = -2RMA~r). 
Summing up all the iterations above, we have 

i=o 

= -i?jv,i(TO-r-)(m + r + l), (Bl) 



which gives Eq. (fT&t . We can verify this result by summing 
up instead of summing down, i.e., with the condition a}. = 
and the iterations 

i^-^o^ilif = 2i?A.,(™+l) 

{ai^-l^f ~ {a^;^f = 2RN,ir 
we have, after adding them up, 

r — m 

= 2i?^,,^(TO + j) 
= Rn ,i{r — m + l){m + r) , 
which is the same as Eq. (IBlb . 



Appendix C: Deriving the excitation spectrum of the 
superconducting qubit chain 

1. State vector compositions for a general W-qubit 
superconducting qubit chain 

The form into which the system Hamiltonian is split as in 
Eqs. (fTTIi-lfTSll ensures that [H^^Hi] = 0. The commutation 
of these two parts implies that we can find simultaneous eigen- 
vectors for iJo and Hi . 

First, for an eigenvector |n;r, to) (or written as \u,r)) of 
-ffo, we have 

Hq |n; r, m) = Hq \u, r) = ujq{m + n) = ujqU 

where {u, n, to} assumes meanings as described in Sec. |IV] 
Note that the eigenvalue ujqU is degenerate, for different com- 
binations of m and n that add up to the same u. Therefore the 
eigenstate of i/o can be written as a superposition 

|u,r) = Cn \n] r, m) S{u — n — to) 

= ^ c„ |n; r, It - 71) A (CI) 

n 

where A is a range delta function 




1, — r < u — n < r 
0, otlierwise 



since we have to ensure the state vectors satisfy the addition 
rules of angular momentum. 

Our next step is to find those of Eq. dClb that are also si- 
multaneous eigenvectors of Hi . With the modified operation 
rule Eq. (fTSI l and setting m = u — n, we can apply Hi to the 
expression and, after reshuffling the terms in the summation 
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such that vectors with the same total excitation number are 
grouped together, we find 

Hi \u, r) \ CnUJon + Cn+nW [n + l)ai'^2„_i 

n 

+ c„_i?7v/na^''2„| |n;r,w - n) A. (C2) 

Since ly is the eigenvalue of Hi, we have 

Hi \u, r) = V |u, r) = CnV |n; r,u — n) A (C3) 



Then comparing Eq. iC2i with Eq. (IC3l l. we deduce a differ- 
ence equation 

c„+i'n\/ (n + l)ai''2(„+i) - c„5„ + c„_i?7Vna[[i„ = 
where v„ = (v — ujon)/ri and the initial conditions are 
c-i = 0, 

Cu+r+l = 0. 

In addition, from the definition of the A function, n < u + r. 
Now write 



(C4) 



and we have a simplified difference equation 

Cn+l - CnVn + Cn-lU (a^^-n) = 0- 

To find the solution, we multiply each equation starting with 



CtiVfi Cn— 

iw„u„_i + {n- l)u„C„_2 



71—1 1 2 

C2 n - n + ^0 n = 

i=2 J = l J=2 

n— 1 n— 1 

n - ^0 n = 

i=i j=o 

Then with the terminating conditions Ci = CqVo and Co = 1, 
we can sum up the equations to eliminate the middle terms and 
obtain 

n-2 

Cn - Qo^n-l + l)CjaLo + l)<3j+2,«-l = 

where we use a shorthand notation 



Qo,n-l ~ Y\. 



To find the analytical expression for C„, we recursively ex- 
pand the factor Cj 



n-1 
n-2 

= Qo^n-l - X! Q0j-lQj+2,«-l(j + + 

i=o 

n-2 J-2 

+ Qj+2.n-lQo,fc-l(3fc+2,i-l 
j=0 fc=0 

By observing that (3oj-iQj+2,7i-i = Qo,«-i/«jW_,+i and so 
on for each pair of Q's in the terms of each recursive expan- 
sion, we can recursively factorize out Qo.n-i and arrive at 

L"/2J 

J'=o 0i---jfc---iL™/2j> 



1) 



a. 



(r) 



where {ji ■ ■ ■ jk ■ ■ ■ jin/2\) is the index set described in 
Sec. HVl Finally, substituting the above expression back to 
the transformation Eq. (IC4l l. we can obtain the coefficients of 
the excitation eigenvector as in Eq. (|2TI) . 



2. One-excitation spectrum for a 4-qubit superconducting 
qubit chain 

The state vector for the one-excitation 4-qubit SQC (u~\, 
= 2, and m E {—2, —1, 0, 1}) can be written as 



r = 

|u,r) =co|0;2,l 



ci|l;2,0)+C2|2;2,-l) + C3|3;2,-2). 

If we assume cq = 1 as a common factor, the rest three coef- 
ficients can be written as 

ci = i?-i/Vi,o 

C2 = i? ^^2,0 ~ '^2,1 

C3 = R ^^^^S.O — R ^^^"^3,1- 

After plugging in the expression according to Eq. ( l22b . we 
obtain the expressions shown in Table. H] 

To find V, and hence the excitation energy, consider the dif- 
ference equations 

C4-C3V3 + 3C2(aL1)' = 
C3-C2V2 + 2Ci(aL'J)' = 
C2-CiSi+Co(4'^)2 = 0. 

■om the last equa- 



Since C4 = and Ci = vqCq, we derive fr 
tion 



C2 = Co koh - {a^o^f 
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and from the first equation 

C3 = 3Co 



and tlie solutions are 



V3 



Substitute these expressions into the second equation and with 
Eq. ( fT6b , we have 

12(z;oWi + 6Rn,i)Rn.i - {vovi - 6Rn,i)v2V3 

+ URNjVoi's = 0. 

Expanding the v, we arrive at a fourth-order polynomial equa- 
tion 

- 6CjoV^ + [llwi^ - SORnjv'^] 

- 6 [wg - 13a;oi?Ar^;?7^] v 

- 36RnjV^ [wq + ^Rn.iV^] = 0. 

If we consider the case with loq = (the conventional TC- 
model case), we have 

~ SORnjV^v^ - 72R% ir]'^ = 

and the roots are 



V = ±yj {15 + 3V33)Rn,iV 

On the other hand, if we consider a weak coupling case ujq ^ 
then the equation becomes 



f"' — Gwqw'^ + IIuiqV^ — 6ijJqV — SGRNjifi^o = 



,2,~,2 



Hence, the excitation energy can be written as in Eq. ( |23l ). 

Now the coefficient Ci for the highest eigenenergy state is, 
since usually ujq < 0, 



^/QR^m 2 2V 



where VL = ujq/ [rj^JORiqj). Since Rn.i < 1, 



\n\ > 



which means that the deformation leads to a larger coefficient 
Ci than that of the undeformed case. For C2, we have 



C2 



since cjq ^ 77 and so |ri| ^ 2. This means that C2 increases by 
a greater proportion than ci due to the deformation. Similarly 
C3 increasing by an even larger factor. The change in the co- 
efficients shows that the a larger deformation favors the states 
with larger number of photons and a more ordered spin-chain. 
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